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Distribution of the eigenvalues of a random system of 
homogeneous polynomials 


Paul Breiding ^ Peter Biirgisser ^ 

Abstract 

Let / = (/i, be a system of n complex homogeneous polyno¬ 

mials in n variables of degree d. We call A G C an eigenvalue of / if there 
exists V G C"\ {0} with f(v) = Xv, generalizing the case of eigenvalues of 
matrices (d = 1). We derive the distribution of A when the fi are indepen¬ 
dently chosen at random according to the unitary invariant Weyl distribution 
and determine the limit distribution for n —^ oo. 

AMS subject classifications: 15A18, 15A69; 60D05 

Key words: tensors, eigenvalues, eigenvalue distribution, random polynomials, 
computational algebraic geometry 

1 Introduction 

The theory of eigenvalues and eigenvectors of matrices is a well-studied subject 
in mathematics with a wide range of application. However, attempts to generalize 
this concept to homogeneous polynomial systems of higher degree have only been 
made very recently, motivated by tensor analysis [12, 13], spectral hypergraph the¬ 
ory [9] or optimization [10]. An overview on recent publications can be found in 
[11], where the authors use the term "spectral theory of tensors". 

Following Cartwright and Sturmfels, who in [4] adapt Qi’s definition of E- 
eigenvalues, we say that a pair (n, A) G (C"^\ {0}) x C is an eigenpair of a system 
/ := (/i,..., fn) of n complex homogeneous polynomials of degree d in the vari¬ 
ables Xi,..., Xn if f{v) = Xv. We call v an eigenvector and A an eigenvalue of 
/. If in addition v'^v = 1, we call the pair [v, A) normalized. 

By [7, Theorem 1.3] we expect the task of computing eigenvalues of a given 
system to be hard. It is therefore natural to ask for the distribution of the eigenval¬ 
ues, when the system / is random. 
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In the case d = 1 we obtain the definition of eigenpairs of matrices. In [5] 
Ginibre assumes the entries of a complex matrix A = {aij) to be indepen¬ 
dently distributed with density exp(—|ajj p) and describes the distribution of 
an eigenvalue A, that is chosen uniformly at random from the n eigenvalues of A. 

Can Ginibre’s results be extended to arbitrary degree d? The answer is yes and 
provided in this paper. Let us call two eigenpairs {v, A), {w, if) equivalent if there 
exists some t € C\ {0}, such that {v, A) = {tw, t'^~^rf. Note that if both (n, A) 
and {w,rj) £ C are normalized, then we must have |f| = 1. This implies that the 
intersection of an equivalence class with the set of normalized eigenpairs of / is a 
circle, that we assume to have volume 27r. Cartwright and Sturmfels point out in 
[4, Theorem 1] that if d > 1, the number of equivalence classes of eigenpairs of a 
generic / is D{n, d) := {dP — l)/(d — 1). 

We define a probability distribution on the space of eigenvalues as follows: 

1. For each 1 < i < n choose fi independently at random with the density 

TT~^ exp(— ||/i|p), where k := Here || || is the unitary invari¬ 

ant norm on the space of homogeneous polynomials of degree d defined in 
Section 3, see also [3, sec. 16.1]. (The resulting distribution of the fi is 
sometimes called the Weyl distribution.) 

2. Among the D{n, d) many equivalence classes C of eigenpairs of /, choose 
one uniformly at random. 

3. Choose an normalized eigenpair {v, A) G C uniformly at random. 

4. Apply the projection (n, A) i—)• A. 

We denote by : C K>o, A i—)• p”’'^(A) the density of the resulting 
probability distribution. Observe that if d = 1, then the density of Ginibre’s 
distribution. 

The unitary invariance of || |p implies that only depends on |A|, but not 

on the argument of A. We therefore introduce the following notation: 

i?:=2|A|^ (1.1) 

We will prove that the random variable R follows a distribution that, if d = 1, is 
mixed from x^-distributions with weights from the uniform distribution on n items, 
and, if d > 1, is mixed from x^-distributions with weights from the geometric 
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distribution Geo{p) truncated at n. (See (2.1) for details on the truncated geometric 
distribution.) 

Here is our main result. 

Theorem 1.1 Let n,d > 1 and A be distributed with density Let denote 
the density of R = 2\\\‘^. 

1. If d = 1, then 


*■(«) = -Edi(«)=E, = 


k=l 


k=l' 


2. > 1, then 




k=l 


V Prob {X = k\ X <n} xlkiR)- 

^X~Geo(l-i) 


Here X 2 k(^) ;= (e 2 ^)/{2^{k — 1)!) is the the density of a chi-square dis¬ 

tributed random variable with 2k degrees of freedom. 

We note that Prob |Af = A;} is the probability that the first success of in- 

X~Geo(p) 

dependent Bernoulli trials, each with success probability p, is achieved in the A;-th 
trial. Moreover, for 1 < A; < n and 0 < (? < 1 we have 

00 

t=0 

Prob {X = kl X <n}; 

X^Geo(l-q) 

for the last equality see (2.1). One can therefore sample |Ap by the following 
procedure. 

1. If d = 1, choose A: G {1,..., n} uniformly at random. 

2. If d > 1, make Bernoulli trials with success probability 1 — ^ until the first 
success. Let £ be the number of the last trial and k the remainder of £ when 
divided by n. 


Prob {X = k + tn\ = 

X~Geo(l-q) 
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3. Choose xi,..., X 2 k ~ ^(0,1). 

2k 

4. Put R xj. 

i=l 

5. Output: ^R. 


Remark 1.2 By de I’Hopital’s rule we have lim This implies that 

lim which yields a connection between the cases d = 1 and 

d —^ 1 

d > 1 (observe that here we allowed d to be any real number). 


We can compute the expectation of the random variable |Ap; cf. Figure 1.1. 


Corollary 1.3 If d = 1, then E |Ap = Ifd > 1, then 


2 ^ n- (n + l)d + d^+^ 

' (d--l)(d-l) ■ 

We have lim E |Ap = 1 and lim E |Ap = ifd > 1. Moreover, for 
fixed n, the function d i—)• E | Ap w strictly decreasing. For fixed d, the function 
n I—?• E \X\‘^ is strictly increasing. 
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Figure 1.1: The left picture shows plots of d i—> E |Ap for n G {2, 3,5,10}. On the 

Xr^P^A 

right are plots of n I—>■ E |Ap for d S {1, 2, 3, 5}. 

X~p^A 

In order to investigate for large n, we can normalize |Ap by dividing it by its 
expectation. We will, however, divide |Ap by n. While for large n this does not 
make a big difference, the formulas appearing are easier to understand. We will 
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also normalize in the case d > 1. So if d = 1, we put 


n 2n ’ 

and if d > 1, we put 

r ■= 1^1' = R(d-l) 

2 lim E |A|2 4d ■ 

n—^OO 

Making a change of variables from i? to r yields the normalized density, denoted 
by pnorm- In [5] Ginibre notes that in the case d = 1 we have 


plormi'T) = 1 [ 0 , 1 ]( t ') : = 

r)—L ’ J 


1, ifO<r<l 
0, else 


(1.3) 


This means that the distribution of the normalized eigenvalue A/yTr converges 
towards the uniform distribution on the unit ball {xGC||x|<1}. 

Our third result covers the case d > 1. 


Theorem 1.4 Let d > 1 be fixed. For any r > 0 we have 


lim 

n^oo 


Pnofmi^) = 2e 


-2t 


Hence, as n ^ oo, the normalized density pummif) converges towards the expo¬ 
nential distribution with parameter 2 (cfi Figure 1.2). 




Figure 1.2: The picture on the left shows plots of pMmC'r) for n € {10, 50,100} together 
with l[o.i] {r). On the right are plots of pndrm('r) for n G (2,3, 5} together with 2e“^'^. 
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1.1 Relation to prior work 

Our definition of eigenpairs is inspired by the following definition of E-eigenvalues 
of tensors and supermatrices given by Qi in [13, sec. 2-3] and [12, sec. 1]. 

Let <1> be a multilinear map —)• M”, that is represented by the real super¬ 
matrix A = For V G M”, Qi puts := A{v,... ,v) (see [13, sec. 3, 

eq. (7)]) and then defines A G C to be an E-eigenvalue of A if there exists v £ 
such that Av^ = \v and v'^v = 1. This definition of eigenvalue is independent of 
the change of orthonormal coordinates. Therefore, A can be regarded as an eigen¬ 
value of <1> itself. Although assuming <1> over the reals, Qi allows the eigenvalue to 
be complex. If the eigenvalue A is real, he calls it a Z-eigenvalue. 

In [4] Cartwright and Sturmfels relax the definition of Qi by considering order- 
{d + 1) tensors/supermatrices over the complex numbers C. They define a pair 
(u, A) G (C"'\ {0}) X C to be an eigenpair of A, if 

Av’^ = \v. (1.4) 


Observe that Qi’s condition v^v = 1 implies that n / 0, while Sturmfels and 
Cartwright require the eigenpair to be an element in (C\ {0}) x C. In reference to 
Qi they call an eigenpair [v, A) satisfying v'^v = 1 normalized, whereas we call 
an eigenpair {v, A) normalized, if it satisfies v'^v = 1. In fact, A is a system of 
homogeneous polynomials over C in the entries of v. So (1.4) coincides with our 
definition. 

Another approach to define eigenpairs of homogeneous polynomial systems is 
given by Lim [10] in his variational approach, which is as follows. 

We denote by || ||fc the £fc-norm on M" for A; > 1. Suppose that F{X) is a 
real homogeneous polynomial in n variables X = (Xi,..., X„) of degree d -|- 1. 
In order to optimize F on the £fc-sphere {||a;||fc = 1}, one can consider the Lan- 
grangian of the multilinear Rayleigh quotient F(X)/||X||^^^, that is L{X, A) := 
F{X) — (d + 1)“^A (||X||^“'"^ — 1), where A is an auxiliary variable. Then the 
equation VL = 0 gives 


VF(X) = A 


/sgn(Xi)^Xf-i\ 


Vsgn(X0"X^7 


IXIL = 1. 


(1.5) 


Note that VF{X) is a system of homogeneous polynomials of degree d. If the 
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pair {v, A) G {H't’IU = 1} x M is a solution of equation (1.5), Lim calls v an C^- 
eigenvector and A an £^-eigenvalue of the system VF. In particular, if A: = 2, the 
£2 -eigenvalues (n, A) G {||n ||2 = 1} x M satisfy 

VF{v) = An, ||n ||2 = 1. 

If we relax the definition of £ 2 -eigenvalues by allowing (n, A) to be complex, the 
pair (n, A) is an eigenpair of the system VF{X) in our sense. 

The organization of the paper is as follows. After some preliminaries presented 
in the next section, we establish in Section 3 the geometric framework for the 
eigenpair problem. Our concepts and notations are close to the ones from [3, sec. 
16]. We define a probabilify disfribufion on {(/, v, A) | f{v) = An}, fhe solufion 
manifold. The pushforward measure of fhis disfribufion wifh respecf fo fhe projec- 
fion onfo fhe space of eigenvalues is precisely Finally, we prove fhe sfafed 
resulfs in Section 4. 

Acknowledgements The basis of Ibis work was laid during fhe program "Algo- 
rifhms and Complexify in Algebraic Geomefry" af fhe Simons Insfifufe for fhe 
Theory of Compufing. We are grafeful for fhe Simons Insfifufe for fhe sfimulaf- 
ing environmenf and fhe financial support. 

2 Preliminaries 

2.1 Differential geometry 

We denote by {x, y) := x'^y the standard hermitian inner product on C". Fur¬ 
thermore, we set ||x|| := (x, x) and S(C"') := {x G | ||x|| = l}. Given 

some X G C’^\{0} we denote by := {y G C” | (x,y) = O} the orthogonal 
complement of x in C”. 

If M is a differentiable manifold and x G M we denote by TxM the tangent 
space of M at x. 

Lemma 2.1 Let v G S(C^). Then T^S(C^) = {a G | sR(a,x) = O} = 
Ty © Riv and this composition is orthogonal with respect to inner product on 
r^S(C"'), that is induced from { ,). 

Proof See [3, Equation (14.11)] and [3, Lemma 14.9]. ■ 
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If M and N are differentiable manifolds and F: M ^ N is differentiable, we 
denote by DF{x) : TxM —)• Tp(^x)^ its derivative at x € M and by NJ{F){x) 
its normal Jacobian at x. 

For more details on normal jacobians and the coarea formula we refer to [3, sec. 
17.3], 


Theorem 2.2 (Coarea formula) Suppose that M, N are Riemannian manifolds of 
dimensions m, n, respectively. Let ^ M ^ N be a surjective smooth map. Then 
we have for any function x ■ M that is integrable with respect to the volume 
measure of M that 



X 

mnIW) 


tiw-'to) 


dN. 


Let El , E 2 be finite dimensional complex vector spaces with hermitian inner 
product, such that dime Ei > dime E 2 . Assume that we have a surjective linear 
map (f> : El ^ E 2 (think of as a derivative and Ei, E 2 being tangent spaces). Let 
r((/)) := {(x, f{x)) ^ El X E 2 } be the graph of f. Then r((/)) is a linear space 
and the projections pi : r(0) —)■ Ei and p 2 ■ L{(j)) —)• E 2 are linear maps. 

The following result is Lemma 3 in [2, sec. 13.2], combined with the comment 
in Theorem 5 in [2, sec. 13.2]. 


Lemma 2.3 Let W be the orthogonal complement o/kerp 2 - Then we have 


|det(pi)| 

|det(p2|ty)| 


\det{4>(p*)\ ^ 


2.2 Expectation of the truncated geometric distribution 


The geometric distribution with parameter p truncated at n > 1 is defined fo be 
fhe disfribufion of a geomefrically disfribufed random variable X wifh paramefer p 
under fhe condifion fhaf X < n. Ifs densify is 


Prob {X = k\X <n} 

X~Geo(p) 


Prob {X = k} 

X~Geo(p) 


Prob {X <n\ 

X^Geo{p) 


1 - g" 


(2.1) 


where q := 1 — p and A: G {1,..., n}. 








Lemma 2.4 Let n > 1 and 0 < g < 1. Then 
E [X\X <n] = — 

X'^Qeo{l—q) 


- (n + l)q^ + 1 


Proof We have Prob {X = k \ X <n} = ^ k G {1, • 


This implies 


X'^Geo(l—g) 

E [X\X <n] = 

X~Geo(l-g) I - q'^ 


k=l 


. ,n}. 


Observe, that Ylk=i ^ derivative of ^ XX ^1-^ = 9 that 


l-Z"+i 


d /I - Z"+i\ _ nZ”+i - (n +1)Z^ + 1 

^ V 1-Z ) ~ (1 - Z)2 


Henee the elaim. 


2.3 The expected characteristic polynomial of a random matrix 

We say that a random variable z on C is standard normal distributed if both real 
and imaginary part of z are i.i.d eentered normal distributed random variables with 
varianee = ^. The eorresponding density is 

t{z) ■■= - exp (-|2:p) , 

TT 

and we write z ~ N{0,^) for this distribution. The reason why we have put 
^ is that for a gaussian random variable z ~ iV(0, on C we have 

E = 1. (2.2) 

Suppose that i? is a finite dimensional eomplex veetor spaee with hermitian 
inner produet and let k := dime E. We define fhe sfandard normal densify on fhe 
spaee E as 

Te{z) ■■= ^ exp (-||z|p) . (2.3) 

if is elear from fhe eonfexf whieh spaee is meanf, we omif fhe subseripf E. Lef In 
be fhe n x n identify mafrix. If a eomplex mafrix A G is disfribufed wifh 

densify (/?cnxn, we wrife A N{0, ^E)- 

Reeall fhaf fhe Gamma funefion is defined by r(n) := for a 
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positive real number n > 0. It is well known that r(n) = (n — 1)! if n is a positive 
integer. The upper incomplete Gamma function is defined as 

roo 

T{n,x) := / 

J t=X 

where x > 0. 

Lemma 2.5 

n—1 ^ 

1. Let X > 0 and n > 1. Then r(n, x) = (n — 1)! e~^ 

fc =0 

2. We have E |det(A)p = n! = r(n — 1). 

3. For I C [n] ;= {1,..., n} we define Aj G to be the submatrix of 

A G indexed by I. Then for any t € C we have that det {A + tin) = 

E/c[n] det 21/. 

Proof The first assertion is from [6, p. 949], the second is [3, Lemma 4.12], and 
the third assertion is a well known fact, cf. [8, Theorem 1.2.12]. ■ 

Proposition 2.6 We have for A G and f G C 

E |det(74 + f/n)P = (n + 1, |fp) . 

A~Ar(o,|r„) 

Proof By Lemma 2.5(3), det(^ + tinxn) = Z^oefo i}” det Aa, hence 
|det(^ + f/„xn)|^ = det^do, detAp. 

Due to Lemma 2.5(2), we have E [det Aa det Ajfi = |q!| !, since we deal with 

centered distributions. Hence, 

E |det(^ + f/nxn)P = ^ k\ r(n + 1, |fp); 

k=0 ^ ' 

the last equality by Lemma 2.5(1). ■ 
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3 Geometric framework 


3.1 Eigenpairs of homogeneous polynomial systems 

Let n,d > 1. We denote by := TidiXi,..., Xn) the vector space of homo¬ 
geneous polynomials of degree d in the variables Xi, ..., Xn over the complex 
numbers C of degree d. The Bombieri-Weyl basis is given by the := (^) 

\a\ = d. We define an inner product on T-Ld via 

^ ^ ^ ^ ^ ^ ^ Claba' (3-1) 

a 06 cx 

The product (3.1) extends to {'Hd)^ in the following way. Let / = (/i,..., /„) 
and £/ = (pi,..., gn) G {UdT- Then we define (/, g) := YA=i{fi^9i)- Moreover, 
for / G {UdT we sef ||/|| := V(/, /). 

Remark 3.1 1. The inner producf (3.1) is fhe unique unifary invarianf producf 

on T-Ld (up to scaling). See [3, Theorem 16.3] and [3, Remark 16.4]. 

2. Suppose that / = (/i, ...,/„) e {UdT and /* = Xla ^i^aea, I < i < n. 
Let k := dimUd and put A := (oj^Q,) G Then ||/|| = ||^||f. where 

II ||f is the Frobenius norm. 

For the sake of clarity, we recall the definition of eigenpairs given in the intro¬ 
duction. 

Definition 3.2 An eigenpair of / G (T-Ld)^ is a pair (v, A) G (C"'\{0}) x C such 
that f{v) = Xv. We call v an eigenvector and A an eigenvalue of /. Further, we call 
eigenpairs {v, A) and {w, rj) equivalent, (u, A) ~ {w, g), if there exists a nonzero 
t G C such that [tv, = {w, g). 

We already noted that the number of equivalence classes of a generic system / 
equals D{n, d) = {dA — l)/(d — 1) if d > 1, cf. [4]. 

3.2 The solution manifold 

Let A := C[Xi,..., Xn, A] be the space of polynomials in the n + 1 variables 
Xi,...,Xn,X We consider the map F : {UdT A^,f f{X) - AX. 
For / E {'Hd)"' we set Fj := F{f), such that 

Ff.C^xC^ C", (u, A) ^ fiv) - Xv. (3.2) 
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Observe that Ff{X,A) consists of two parts, one homogeneous of degree d and 
one homogeneous of degree 2. Let us denote by dx and Oa the partial derivatives 
of Ff{X,A) with respect to X = (Xi,... ,X„) and A, respectively. Then the 
derivative of Fj at {v, A) has the following matrix representation: 


dxf — dx{AX), 


-9a(AX) 


iX,A)=iv,\) 


dxf{v) - XIn, -V 


(3.3) 


where denotes the n x n-identity matrix. 

We adapt the terms “solution manifold” and “well-posed” from [3, sec. 16.2] 
and tailor them to our (structured) set {Ff \ f G {'Hd)"'}- Compare [3, Open Prob¬ 
lem 15]. We call 


V := {(/, n, A) G {ndT X S(C-) X C | Ff{v, A) = 0}, 


the solution manifold and its subset 

W ;= {{f,v,X) G V I vkDFf{v,X) = n} 
the manifold of well-posed triples. 

The group ff(n) of unitary linear transformations C" —)• C” acts on {T-Ld)^ and 
V, respectively, via 

U.f ■.= UofoU-^ and U.{f,v,X):={U.f,Uv,X). (3.4) 

We note that W is invariant under the group action and that U (n) acts by isometries; 
see [3, Theorem 16.3]. 

Lemma 3.3 The solution manifold V is a connected and smooth submanifold of 
{'Hd)"' X S(C"’) X C of dimension dimjRV = dim]R(?fd)’^ -|- 1. Moreover, the 
tangent space ofV at (/, v, A) equals 

{(/,(;, A) G fHdT X T,S(C") x C | f{v) + DFf{v,X){v,X) = o} . 

Proof The map G : {'Hd)"' x S(C”) x C —)> (/, v, A) i-)- Ff{v, A) has V as 

its fiber over 0. The derivative of G, 

DGif,v,X): (?fd)"xr,S(C")xC^C", if,v,X)r^fiv)+DFfiv,X){v,X), 
is clearly surjective. Therefore 0 G is a regular value of G and Theorem A.9 in 
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[3] implies the assertion. 


The following lemma is easily verfied using Euler’s identity for homogeneous 
functions. 

Lemma 3.4 Let {f,v,X) G W. Then ker DFf{v, X) = C{v,{d— 1)A)^. In 
particular, DFf{v, A)|r„xC invertible. 

Corollary 3.5 The tangent space T(^f^y^x)V at (/, A) G W is given by 

[{f,v + riv,X) G {TLdT x {T^®Riv) x C | ('u,A) = -DFf{v,X)\:^^^^f{v)y 

Proof Let (/, A) G V be fixed. By Lemma 3.3 the tangent space of V at 
(/, V, A) equals 

[{f,w, A) G {TidT X T,E{C^) X C I DFfiv, A) {w, A) = -f{v)] . 

From Lemma 2.1 we know that r„S(C”) = © Mir;. Lemma 3.4 tells us that 

keiDFf{v, A) = C(r;, (d — 1)A) = M.{v, (d — 1)A) © Mi{v, {d — 1)A). 

Hence, {TyE{C^) x C) n ker DFf{v, A) = {d — 1)A). From Lemma 3.4 we 
know that DFf{v, A)|7’„xC is invertible. We conclude that if (/, w, A) G T(j 
then there exist uniquely determined i) ^ Ty and r G M such that w = i) + irv and 
DFf{v, X){w, A) = DFf{v, X){v, A), from which the claim follows. ■ 


3.3 Projections and normal jacobians 

We consider the projections 

vn: V^(Hdr,(/,r;,A)^/, vrs : V ^ S(C") x C, (/, r;, A) ^ (r;, A). (3.5) 


It is essential that the quotient of the normal jacobians of tti and tt 2 can be 
computed in the following way. 


Lemma 3.6 For all (/, u, A) G W we have 


NJ{Tri){f,v,X) 

NJiTr2)if,v,X) 


\det{DFf{f,v,X)\Tyxc)f ■ 
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Proof Let U{n) be group of unitary maps C” C". Recall from (3.4) that 
for U G U{n) and {f,v,X) G V we have put U.{f,v,X) := {U.f,Uv,X). By 
definition, the projections vri, tt 2 are (n)-equivariant. Hence for any U G U{n) 
we have N J{7ri){f,v, X) = N J{TTi){U.{f,v, X)), i = 1,2. It therefore suffices to 
show the claim for v = ei, where ei := (1, 0,..., 0) G C”. 

Suppose that (/, ei, A) G W. The derivatives of tti and tt 2 are the projections 

D7nif,ei,X): {f,v,X)^f, 

D7r2{f,ei,X): ^ re,§(C") x C, {f ,v, X) ^ (v, X). 

Let us write f = where for all a we have fa G C^. Then we obtain 

a 

/(ei) = f(d,o,...,o)- Hence, / /(ei) is an orthogonal projection. We will denote 

this projection by H. By Lemma 2.1 the projection i) + riv i—)• u is orthogonal 
as well. Using 3.5 it follows that is the graph of the surjective linear 

function 


-DF/(ei, ° (n X 0) : {HdT x Te^ x C. 

Applying Lemma 2.3 yields the claim. ■ 

3.4 The eigendiscriminant variety 

We define the set of ill-posed triples (/, v, X) to be 

S' := {(/,u,A) G V I ikDFf{v,X) < n} = V\W. (3.6) 

Moreover, we denote by S the Zariski closure of 7ri(S'). In reference to [1], we 
call S the eigendiscriminant variety. 

Remark 3.7 We have (/, v, A) G S' if and only if (u, A) is not an isolated root of 
the polynomial Ff. Thus, / G 7ri(S') if and only if Fj has a double root or / has 
infinitely many roots. 

Proposition 3.8 1. We have / 0 S, if and only if the number of equivalence 

classes of f equals D{n, d). 

2. The set S is a closed hypersutface of (Hd)^ of degree at most n(n — 
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Proof For Item 1 use [4, Theorem 1.2]. In [1, Theorem 4.1, Corollary 4.2] it is 
shown that the eigendiscriminant variety for tensors in is an irreducible hy¬ 

persurface. We obtain S by intersecting this with the linear subspace of symmetric 
tensors and requiring ||n|| = 1. The assertion Item 2 follows from the dimension 
theorem, Bezout’s theorem and the fact that E is properly contained in {'Hd)^ (see 
3.9 below). ■ 


In [14] the following explicit element in {'Hd)‘^\Tj is described {d > 1). 
Proposition 3.9 Let (l){X) := {Xf, X^ ...,X^)g {UdT. Then 


Proof One has 


F^iX,X) 


/xf-AXi\ 

\X^-Axj 


We are going to show that (p has exactly D{n, d) many classes of eigenpairs. 
Clearly, for any v G C\{0} we have F^(u,0) / 0. Hence, any equivalence 
class of eigenpairs of (p contains some representative of the form (u, 1). Let be a 
primitive {d — l)-th root of unity and define 


M := {(eiCS ..., enC") | e C {0, 1}^ \ {0} , Vj : 1 < z,- < d - l} 


Observe that F^{z, 1) = {zf — Zi)2^i = 0, if and only if z G M U {0}. For 
all z G M and t G C we have {z, 1) ~ {tz, 1) if and only if f = C* for some 
1 < z < d — 1. Let if := {() denote the cyclic group generated by We 
define a group acfion of if on M via componenfwise mulfiplicafion. The number 
of equivalence classes of eigenpairs of (p fhen equals fhe number of il-orbifs in M. 
For zz G if puf M“ := {z £ M \ uz = z}. Observe fhaf for zz / 1 we have fhaf 
M“ = 0. Using Burnside’s lemma we obfain 

number of il-orbifs in M = ■;—r = ■;—r \M\ = —-= Din, d). 

if ^ n d-l 

' ' «eit ' ' ■ 


3.5 The standard distribution on the solution manifold 

The definition of standard distribution is adapted from [3, eq. (17.19)]. Following 
(2.3), we say that a random variable / on {TLd)'^ is standard normal distributed, if 
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/ has the density 


‘fif) ■= ^iHd)"if) = ^ exp (-||/f) , where k = dimc(?^d)"- 

By construction (p{f) is invariant under the action of Z^(n). 

The following procedure: 

1. choose / according to the standard normal distribution. 

2. choose some normalized eigenpair (u, A) of / uniformly at random. 

yields a probability distribution on V, which we call the standard distribution and 
denote it by (/, v, A) ~ STDy. Clearly, the standard distribution is invariant under 
the action of U (n) on V. 

Observe that the two steps above are precisely the steps Item 1-Item 3 in the 
operative description of given in the introduction. This implies that p^’'^ equals 
the density of the pushforward measure of STDy with respect to the projection 
TTs: V C, {f,v,X) ^ A. 

According to 3.8, the fiber 

V{f) := {(u, A) G S(C") X C I if,v,X) G V} = 7 r2{7T^\f)) 

over / 0 S consists of iA = D{n, d) disjoint circles, each of them having volume 
27r. Hence the density of the uniform distribution on V (/) equals (27rD)~^. As in 
[3, Lemma 17.18], one can now show that the density of the standard distribution 
is given by 

PSTDv(f,v,X) = NJ{7ri){f,v,X)ip{f), (3.7) 

where tti : V —)• (Tfd)” is the projection from (3.5). 

We denote by 

V(v, A) := {/ G (Tidr I (f, V, A) G V} = 7ri(7r2'^(u, A)) 
the fiber of 7r2 over (v, A) G S(C”) x C. 

Lemma 3.10 Let 0 : V —)■ M. be an integrable map that is invariant under the 
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group action from (3.4) and ei := (1,0,..., 0) G S(C"). Then 




V, A) psTDv if, V, A)dV = 


TT 


71—1 


T{n)D{n,d) 


E{X)dC. 


AeC 


where 

E{\)= I \det DFf {ei, X)\^ 9{f,ei,X)pvie„x}if)dViei,X). 

/6V(ei,A) 

Proof Using the coarea formula, we obtain 


i 

'(D,A)e§(C")xC 


A)pstDv(/,^^, A)dV 


f <9(/,7;, A)/3stDv(/,^, A) 

lf£V{v,\) NJi7r2)if,v,X) 


dV{v,X) 


d(S(C") X C) 


By the definition of /9 stDv> Lemma 3.6, and the unitary invariance of 9 we have 
that 


&{f,v,X) 


lfeV{v,x) NJ{7r2)if,v,X) 


PSTDv dU (v, X) 


2 TrD 


2TrD . 

EjX) 

2ttD' 


'f&V{v,X) 


|det DFf{v, A)|t„xc|^ dif, v, A) ip{f) dV (n, A) 


UeV'(ei,A) 


|det DFf{ei, X)\t,^ xc|^ dif, ei, A) ipif) dU(ei, A) 


(3.8) 


Observe that the integral (3.8) does not depend on v anymore. The claim follows 
by using f ldS(C") = ^ . 


4 Proofs 

We are now ready to prove 1.1. 

Proposition 4.1 The pushforward density o/STBy with respect to is 

= 'r ‘-'r (n. If) ^ ]_( 

7rD{n,d) L(n) TrD{n,d) \ d J 
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Proof Before we start, we remark that Lemma 2.5 justifies the right equality. By 
Lemma 3.10, the pushforward distribution is obtained by computing 


y./ ■.p.f n / detDF/(ei,A)|T,jxC ^v{eu\){f)dV{ei, X) (4.1) 

r(n)D(n,d) 7jgv'(ei,A) 

The case n = 1 is an easy exercise. So let us assume that n > 1. Observe that 
V (ei, A) is the affine space 

y(ei,A) = AXfei + {g £ {UdT I 5(ei) = o}. 

Let 22 := |/i G {'Hd)"' \ h{ei) = 0,Dh{ei) = o|. By [3, equation (16.10)], for 
/ £ ^(ei, A), there exist uniquely determined h ^ R and M G 
that we can orthogonally decompose / as 

/ = XXfei + X^-^VdMX' + h, (4.2) 


where X' = (X 2 , ■ ■ ■, Xn)'^. We have that 


dxf{ei,X) 


dxJ{ei,X), dx'f{ei,X) 


dXei, 


VdM 


G C" 


(4.3) 


Let a G be the first row of M and A G be the matrix that 

is obtained by removing the first row of M. By (3.3) and (4.3) the derivative of Fj 
at (ei, A) has the matrix representation 


dxf{ei,X)-XIn, -ei 


{d — 1)A Vda 
0 VdA — XIn-i 


-1 

0 


g (j-<nx(n+l) 


This implies det iAFy(ei, Ajl-r^j xC = — det {y/dA — XR-i)■ 

The summands in (4.2) are pairwise orthogonal. From this we get that ||/|p = 
|Ap + ||M|||, + ||/i|p, which implies that 




=-|Ap 


TT" 
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Integrating over a and /i in (4.1) therefore yields 


Jf&V{eu\) 

=-|Ap 


I det DFf (ei , A) I Tej X c I V^v'(ei ,A) (/)df^(ei , A) 


E 

(^n-lg-lAp 


det (^VdA — A/„_i^ 


TT 

Tj-n 


E 


det ( A - 




the last line by 2.6. Plugging this into (4.1) the claim follows. ■ 

Proof (Proof of 1.1) 4.1 shows that the distribution of the eigenvalue A only 
depends on |A|. As in (1.1) we put r := |A| and R := 2r^. Making a change of 


variables, we obtain the density p^'^{R) := | p^’^{r). From 4.1 we obtain 


_TT ^nA{ 


r! 

nA/ j^\ ^ 


n—1 


2D{n, d) 


^1-1 / D 'v ^ 

e-f 

^ k\ \2d 


k=0 


(4.4) 


If d = 1, (4.4) becomes 




e-2R>^ 


n 

-E 

71 ^ 


g 2 1 


n ^ 2^+iA:! n ^ 2^{k - 1)!' 
fc=0 k=l ^ ’ 

For any k we have that e~^ R^~^/{2^{k — 1)!) is the density of a chi-square 
distributed random variable with 2k degrees of freedom, which proves the assertion 
in this case. 

If d > 1, put q := 2, such that D{n,d) = (1 — — q))- Then (4.4) 

becomes 


6 ^ (1 — q)q^ 

* ~ 1-9” ~ A 2‘’{fc - 1)! 1 - 9' 


E 


e 2 1 (1 — q)q 


k-l 


Using that Prob {X = k \ X < n} = q^ ^(1 — q)/{l — q^), see (2.1), fin- 

X~Geo(l—g) 

ishes the proof. ■ 


To prove 1.3 we will need the following lemma. 

Lemma 4.2 Let n > 1. 
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1. Ifd = l,then E lAP = E [X], 

2. Ifd > 1, then E |A|2 = E [X\X <n]. 

\^p^,d X~Geo(l-i) 

Proof We prove the claim for d > \. (The case d = 1 is proven similarly.) If 
i? = 2|Ap, then E [|Ap] = i E[i?]. From 1.1 we get 


poo 

E[i?] = / Rp^'^{R)dR 
Jr=o 

^ pn 

= V Prob {X = k\X<n} / RxlkiR)dR. 
^X~Geo(l-i) Jr=o 


= V Prob {X = k\ X <n} 2k = 2 E [X\X<n], 

^X~Geo(l-i) x~Geo(l-i) 


where we have used that a "distributed random variable with 2 k degrees of 
freedom has the expectation 2 k. ■ 


Proof (Proof OF 1.3) Ifd = 1, from Lemma 4.2 we immediately get E |Ap 

n+1 

2 • 

If d > 1, by Lemma 4.2, we have that E |Ap = E [X \ X < n]. 

Xr^p^.d jV~Geo(l-i) 

Therefore, Lemma 2.4 with q '■= \ implies 

2 _ n - (n + l)d + d’"+^ 

A~p",^' ' (d--l)(d-l) 

as claimed. For fixed n we obtain 

, n — (n + l)d + d”+^ n + 1 

hm ——-r— -;— = - 

d^i (d"-l)(d-l) 2 

by using de I’Hopital’s rule twice. Therefore, the map 


M>i —>• M, d I—>■ 


71+1 
2 ’ 

n-(n+l)d+(i"+l 


if d = 1 
if d > 1 


is continous and differentiable on M>i. One checks that its derivative on ]R>i is 
negative. Hence, for fixed n, we see that d i—)• E [| Ap] is strictly decreasing. In 

the same way we can prove that, if d is fixed, n i—)> E [| Ap] is strictly increasing. 

A~p".d 
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Further, 


, n — (n + l)d + d , ng"'+^ — (n + l)g + 1 

lim ——---;— = lim-;--r-= 1 

d-^oo {d'^ — 1) {d — 1) g-s-o {1 — q^) (1 — q) 

If d > 1, we have 

lim E [|Afl = lim ~ *"iT ‘ . 

n^oo Xr^pn.d n-Aoo (1 — q"^) [1 — q) ^ ~ Q 

where again q = 

Proof (Proof of 1.4) Let d > 1. Recall from (1.2) that we have put r 
and that we denote the density of r by pnorm- Using 1.1 we get 


P^.iLir) = 


Ad 


n,d 


Adr 


\d-lj d^-1 
Again putting g = ^ we obtain 


2d^ 1 / 2r ^ 

e d-i y — 


E 

A:=0 


kl \d-l 




1 - 


k=0 


kl \1 — q 


Since 0 < g < 1, we have lim g” = 0. Hence, 


which finishes the proof. 


oo „ \ k 

^ 1 f 2qT \ _ 2r 

— z e , 


k=0 
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